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Preface 


The  work  presented  In  this  report  is  a  combination  of  the  investigations  pur* 
sued  under  three  different  in-house  Work  Units  being  conducted  by  the  Cloud 
physics  Branch  of  AFGL.  The  measurement  and  determination  of  actual  hydro¬ 
meteor  particle  size  distributions  in  the  free  atmosphere  is  a  difficult  basic 
research  problem1  and#  as  such,  is  being  pursued  under  Task  2310G5,  Cloud 
Physics,  which  is  funded  by  the  Air  Force  Office  of  Scientific  Research  (AFOSR). 
Data  used  to  demonstrate  the  application  of  the  method  were  obtained  while  the 
MC-130E  and  Branch  members  were  supporting  the  SAMSO/ABRES  ANT-2  mission 
at  the  Kwajalein  Missile  Range  under  the  Advanced  Development  Project  627A, 
Weather  Erosion  Program.  However,  the  majority  of  the  effort  reported  herein 
was  in  the  nature  of  exploratory  development  and  was  conducted  under  Task  667012, 
Cloud  Microphysics. 

Thanks  are  expressed  to  Dr.  Arnold  A.  Barnes,  Jr.  and  Mr.  Vernon  G.  Plank 
for  their  valuable  suggestions  in  the  writing  of  this  report. 
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1.  Barnes,  A. A.,  Jr.  (1978)  New  cloud  physics  instrumentation  requirements, 
Preprints.  Fourth  Symposium  on  Meteorological  Observations  and 
^tn^^jttion.  ^AmerlcaiT Meteorological  Society.  284-268  AFQL-Tft-0093, 
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A  Method  to  Predict  the  Parameters 
of  a  Full  Spectral  Distribution 
From  Instrumental  ly  Truncated  Data 


1.  INTRODUCTION 

The  Air  Force  Geophysics  Laboratory  (AFGL)  has  equipped  two  aircraft,  a 
Learjet  and  MC-130E,  with  a  variety  of  instruments  for  measuring  meteorological 
parameters.  One  use  of  the  data  acquired  by  these  aircraft,  either  by  themselves 
or  in  conjunction  with  ground-based  radar,  is  to  determine  the  liquid  water  con¬ 
tent*  of  precipitable  storm  hydrometeors. 

Contained  within  this  specialized  instrumentation  are  two  types  of  Optical 
Array  Probes  manufactured  by  Particle  Measuring  Systems  Inc.  of  Boulder, 
Colorado  which  are  specifically  designed  to  characterize  and  count  the  hydro - 

meteors.  A  one -dimens  tonal  system  (PMS  1-D)  classifies  these  hydrometeors  as 

2 

to  their  physical  size  in  one  dimension  and  also  measures  the  number  density. 

The  two-dimensional  instrument  (PMS  2-D)  provides  an  electronically-produced 

shadowgraph  from  which  a  two-dimensional  size  can  be  measured  and  type  of 

3 

precipitation  particle  can  be  inferred. 

Using  the  data  from  the  1-D,  the  liquid  water  content  (M)  of  the  sampled 
hydrometeors  can  be  determined  by 


(Received  for  publication  7  January  1980) 

^Liquid  water  content  is  defined  in  this  report  as  the  mass  of  water  (rain  or  melted 
ice  and  snow)  per  unit  atmospheric  volume.  Units  are  grams  per  cubic  meter. 

The  references  cited  above  will  not  be  listed  here.  See  References,  page  41. 
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M  =  c 


(1) 


i=n 

£  N.D?gm 
i=l 


-3 


-3 

where  (m  )  is  the  number  concentration  in  a  class  having  D.  (mm)  as  the  mid¬ 
size  diameter  and 

c  =  1  X  10‘3  ^  (2) 

where  pw  (g  cm  )  is  the  density  of  liquid  water.  The  equivalent  of  the  radar 
reflectivity  (Z)  can  also  be  calculated  from  the  same  spectral  data  by 

i=n 

z  =  2  Ni  Df mm6  m’3  • 

i=l 

In  the  case  of  an  ice  hydrometer,  the  D  in  the  equations  is  taken  as  the  equiva¬ 
lent  melted  diameter  of  the  ice  crystal  and  is  defined  by 

D  =  a  (4) 

where  a  and  0  are  values  that  change  in  accord  with  the  crystal  type  and  i  is  the 

4 

physical  dimension  of  the  ice  particle  as  measured  by  the  PMS  equipment. 

With  the  spectral  information  supplied  by  this  instrument  and  the  knowledge 
of  crystal  habit,  the  M  and  Z  of  the  distribution  can  readily  be  determined,  but  a 
problem  develops  when  some  of  the  crystals  grow  to  sizes  that  exceed  the  measur¬ 
ing  capabilities  of  the  instrument  (3.  601  mm  for  the  1-D  system  on  the  Lear  and 
5.402  mm  for  the  system  on  the  MC-130E).  At  this  point,  the  spectral  distribution 
is  instrument  ally  truncated  and  the  derived  M  and  Z  are  valid  only  for  the  portion 
of  the  spectrum  that  was  actually  measured.  The  2-D,  on  the  other  hand,  is  cap¬ 
able  of  recording  parts  of  these  large  particles,  and  gives  some  knowledge  as  to 
physical  size  although  the  number  density  and  size  classification  are  not  so  well 
defined. 

The  purpose  of  this  paper  is  to  present  a  method  of  extrapolating  the  1-D 
truncated  distributions  so  that  a  reasonably  accurate  estimate  of  the  total  M  and  Z 
may  be  found  and  then  describe  how  these  corrected  data  may  be  conveniently 


4.  Cunningham,  R.M.  (1978)  Analysis  of  particle  spectral  data  from  optical  array 
(PMS)  1-D  and  2-D  sensors.  Fourth  Symposium  on  Meteorological  Observa¬ 
tions  and  Instrumentation.  April  10-14,  1978,  Denver,  Colorado,  Preprint 
Vol. ,  pp  345-350. 
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incorporated  in  mathematical  relationships  that  more  closely  represent  the  actual 
hydrometeor  spectra. 


2.  EXTRAPOLATION  OF  A  TRUNCATED  DISTRIBUTION 

The  problems  of  extrapolating  a  distribution  are  twofold:  Which  method  is  to 
be  used  to  extend  the  distribution  beyond  the  upper  truncation  limit  and  what  is  the 
limit  of  the  extension,  or  the  determination  of  maximum  particle  size? 

When  addressing  the  first  problem,  it  is  logical  to  assume  that  the  form  of 
the  number-density  distribution  will  remain  relatively  constant  beyond  the  physical 
size  limitation.  Therefore,  the  trend  of  the  distribution  should  be  established 
from  the  measured  portion  of  the  spectra. 

Figures  1,  2  and  3  show  the  plotted  number  density  vs  the  equivalent  melted 
diameter  for  three  separate  passes  taken  on  July  4,  1978  over  the  Kwajalein  Atoll. 
All  are  measurements  of  hydrometeors  classified  as  large  snow,  one  sampled  by 
the  Lear  and  the  others  by  the  MC-130E,  and  all  three  have  maximum  physical 
sizes  that  are  much  larger  than  the  I-D  instruments  are  able  to  measure.  (In  all 
cases,  the  a  in  Eq.  (4)  is  0.  4  and  the  0  is  0.  782).  These  distributions  can  be  con¬ 
sidered  typical  of  the  experimental  data  collected  by  AFGL  in  past  years  in  that 
they  exhibit  a  systematic  decrease  in  the  number  density  over  the  mid-  to  large- 
size  diameter  range.  The  figures  demonstrate  that  a  straight  line,  on  a  log-linear 
plot,  can  adequately  describe  the  population  for  the  major  part  (approximately 
92  percent)  of  the  total  measured  size  range  of  each  of  the  distributions.  Devia¬ 
tion  from  the  straight  line  occurs  at  or  about  the  junction  of  the  cloud  and  precip¬ 
itation  probes,  with  the  cloud  probe  spectra  following  a  line  having  a  much  steeper 
slope.  The  contributions  of  M  and  Z  from  the  precipitation  probe  data  (Mp  and 
Zp)  constitute  the  majority  of  the  measured  values  and  indicate,  at  least  in  trun¬ 
cated  situations,  that  the  M^,  and  Z^  derived  from  the  cloud  probe  spectra  are 
minor  contributions  as  shown  in  Table  1. 

When  the  data  from  the  precipitation  probe  of  the  1-D  instrument  is  considered 

by  itself,  it  is  apparent  that  the  size  distribution  can  be  reasonably  described  by  a 

function  of  exponential  type  which  agrees  with  the  observations  made  by  numerous 
56789  10 

other  investigators.  '  *  •  '  '  Although  this  occurrence  is  evident  in  the  vast 
majority  of  cases,  an  occasional  sampling  level  will  exhibt  number -density  dis¬ 
tributions  that  do  not  conform  to  an  exponential  fit  and  can  not  be  used  in  the  fol¬ 
lowing  processes. 


Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  41. 
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X  CLOUD  PROBE 
o  PRECIPITATION  PROBE 


Figure  1.  Number -Density  Plot  of  Large  Snow  Hydrometeors 
as  Recorded  by  the  PMS  1-D  on  thtf  Lear  at  an  Altitude  of  7  km 
over  the  Kwajalein  Atoll  on  July  4,  1978  (99  sec  av) 


Figure  2.  Number-Density  Plot  of  Large  Snow  Hydrometeors 
as  Recorded  by  the  PMS  1-D  on  the  MC-130E,  Pass  3  at  an  Al 
titude  of  5. 1  km  over  the  Kwajalein  Atoll  on  July  4,  1978 
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Figure  3.  Number -Density  Plot  of  Large  Snow  Hydrometeors 
as  Recorded  by  the  PMS  1-D  on  the  MC-130E,  Pass  5  at  an 
Altitude  of  6.  6  km  over  the  Kwajalein  Atoll  on  July  4,  1978 
(212  sec  av) 


An  exponential  size  distribution  specifies  that  the  number  concentration  of  the 
hydrometeors  will  decrease  with  increasing  diameter  as 

ND  =  Noe~AD  No‘  m"3  mm_1  (5) 

where  Nq  is  the  number  per  cubic  meter  per  millimeter  bandwidth  and  is  a  con¬ 
stant  of  the  exponential -distribution  function,  A  is  number  per  millimeter  band¬ 
width  and  is  the  slope  of  the  number -density  distribution  and  D  is  the  drop  diam¬ 
eter  in  millimeters  or,  in  the  case  of  ice  hydrometeors,  the  equivalent  melted 
diameter.  The  total  number  of  hydrometeors  in  a  population  between  a  specific 

minimum  (d)  and  a  maximum  diameter  (D.,)  is 

M 

dm 

Nt  =  /  NdD  No.  m-3  (6) 

d 
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Table  1.  Comparison  of  the  M  and  Z  Contributions  Derived 
from  the  Measurements  of  the  Cloud  and  Precipitation 
Probes  of  the  PMS  1  -D 


-3 

M  g  m 

Total 

™  -3 

M  g  m 

Percent  from 

Cloud  Probe 

Mc 

Precip.  Probe 
MP 

Cloud  Probe 

Precip.  Probe 

Lear 

0.0225 

0.3279 

0.3504 

6.4 

93.  6 

MC-130E 
Pass  3 

0.0283 

0.3824 

0,4107 

6.  9 

93.  1 

MC-130E 
Pass  5 

0. 0244 

0.3062 

0.3306 

7.4 

92.  6 

Z  mm 

6 

m 

Total 

Zmm6  m"3 

Percent  from 

Cloud  Probe 
Zc 

Precip.  Probe 

Z 

P 

Cloud  Probe 

Precip.  Probe 

Lear 

0.014 

177.2 

177. 214 

<0.01 

>99.99 

MC-130E 
Pass  3 

0.054 

544.  1 

544.  154 

0.01 

99.99 

MC-130E 
Pass  5 

0.040 

336.4 

336.440 

0.01 

99.99 

„  No  rN  „  -3 

Nt  =  — - -  No.  m 


(7) 


where  r^  is  the  truncation  ratio  of  the  number  of  hydrometeors  contained  within 


the  d  and  limits  to  the  total  N  as 


r 

J  NdD 

,  _ 

N  oo 

J  NdD 


(8) 


r  =  e‘dA  -  e”D^A 
rN  e  e 


(9) 
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The  hydrometeor  liquid  water  content  is  distributed  with  diameter  as 


Md  •  ~  X  10-3  pw  No  e'A  D  g  m'3  mm'1 

and  the  total  M  from  d  to  D,.  is 

M 

dm 

M  =  /  Md  dD  ,  g  m'3 
d 


or 


M  = 


»  x  IQ'3  pw  No  r(4)  rM 
6  A4 


g  m 


(10) 


(11) 


(12) 


where  r(4)  is  the  gamma  function  of  4  and  is  the  truncation  ratio  of  the  liquid 
water  content  contained  within  the  d  and  limits  to  the  total  Ivl  as 


dm 

/  MDdD 

rM  = 

/MD  dD 


(13) 


or 

rM  ■  I  {e'dA[(dA)3  +  3(dA)2  +  6dA+  6] 

-e  M  [(DmA)J  +  3(DM’Ar  +  6DmA  +  6]  }  . 

The  radar -reflectivity  distributed  values  for  the  same  population  are 

v  -  m  n®  6  -3  -1 

Ztj*±  *  N  u  e  mm  m  mm  . 

U  o 


(14) 


(15) 
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.v-Au. ;  .  ■-  *\  »f  .  v‘0;. 


.. . ,  ■  ■ 
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and  the  total  Z  within  the  d  to  D__  limits  is 

M 


Z  =  f  ZD  dD  , 


6  -3 

mm  m 


N  r(7)  r„  • 

Z  O  lu  _ D _ <J 

- - s -  mm  m 

A 


where  r(7)  is  the  gamma  function  of  7  and  r^  is  the  truncation  ratio  of  the  radar 
reflectivity  contained  within  the  d  and  D^j  limits  to  the  total  Z  as 

dm 

/  ZDdD 


r  n  • 

Z  oo 


/  ZDdD 


rZ  =  t! o  {e_dA  t(dA)6  +  6(dA)5  +  30(dA)4  +  120(dA)3  +  360(dA)2  +  720dA+  720] 
-  e  °mA  HDmA)6  +  6<DmA)5  +  30(DmA)4  +  120(DMA)3  +  360(DMA)2 
+  720(DmA)  +  720]  }  .  (19) 

When  Eqs.  (12)  and  (17)  are  solved  for  NQ  they  can  be  equated  as 


_ 6  M  A  _  .  Z  A' 

•  X  I®'3  Pw  r(4)  rM  r<7'  rZ 


and  solved  for  A  as 


6  M  r(7)  rz 
X10-3pw  T(4)r# 
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Substituting  r<4)  =  6,  r(7)  =  720  and  pw  =  1  g  cm"  ,  Eq.  (21)  reduces  to 


A  =  61.2 


(*r 


mm 


-1 


(22) 


Since  both  and  incorporate  A  as  a  term  (Eqs.  (14)  and  (19))  then  a  "trial 
and  error"  method  can  be  used  to  solve  Eq.  (22),  whereby  the  value  of  A  is  ad¬ 
justed  until  both  sides  of  the  equation  are  equal.  Once  A  has  been  found,  Nq  can 
be  determined  using  either  Eq.  (12)  or  (17)  and  the  total  number  of  particles  may 
be  calculated  from  Eq.  (7). 

If  the  M  and  Z  derived  from  the  truncated  precipitation -probe  measurements 
of  the  1-D  instrument  are  used  in  these  equations  with  the  appropriate  diameter 
limits,  where  now  becomes  the  upper -truncation  diameter  D^,,  then  the  dis¬ 
tribution  properties  of  the  spectra  can  be  described  by  a  function  of  exponential 
type.  The  plotted  exponentials  of  Figures  1,  2  and  3  (solid  lines)  were  calculated 
in  this  manner  where  d  =  0.  097  mm  and  DT  =  1.  090  mm  for  the  Lear  and 
d  =  0.  133  mm  and  D^,  =  1.496  mm  for  the  MC-130E  data.  It  is  apparent  from 
these  figures  that  the  solid  line  representing  the  exponential  number  distribution 
establishes  the  trend  of  the  spectrum.  An  estimation  of  that  portion  of  the  popu¬ 
lation  missed  because  of  instrument  truncation  can  be  made  by  extrapolating  the 
line  to  some  larger  diameter. 

The  second  part  of  the  problem  is  the  establishment  of  a  new  value  which 
reflects  the  maximum  size  of  the  nontruncated  spectra.  Past  studies  of  exponen¬ 


tial  distributions  conducted  at  AFGL 


have  shown  that  D-.  can  be  related  to  A  as 
M 


DMA  =  c 


(23) 


when  maximum  particle  sizes  of  known  crystal  habit  are  measured  with  the 
Aluminum  Foil  Hydrometeor  Sampler11  and/or  the  PMS  2-D  and  are  converted  to 
Dm  va*ues  though  Eq.  (4).  The  normal  range  of  C  is  between  9  and  12  and  is 
dependent  on  crystal  density  and  configuration.  For  the  large -snow  category  of 
Figures  1,  2,  and  3,  C  has  the  value  of  9  and  was  determined  from  measurements 
of  the  PMS  2-D  and  the  use  of  0.  4  for  a  and  0.  782  for  0  in  Eq.  (4). 

The  maximum  physical  size  recorded  by  the  MC-130E  2-D  for  the  flight  passes 
plotted  in  Figures  2  and  3  were  approximately  8  and  6.  5  mm.  The  calculated  D^ 
in  Figure  2  is  2.0  mm  which  translates  into  a  I  of  7.  82  mm  using  the  reverse 


11.  Church,  J.F. ,  Poes,  K.  K. ,  and  Spatola,  A.  A.  (1975)  The  Continuous 

Aluminum  Foil  Sampler;  design  operation,  data  analysis  procedures,  and 
operating  instructions.  Instrumentation  Papers.  No.  235.  AFCRL-TR- 
75-0370,  AD  A019  630,  70  pp. 
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form  of  Eq.  (4).  The  of  Figure  3  is  1.  698  mm  and  the  corresponding  i  is 
6.  35  mm.  Unfortunately,  no  comparison  can  be  made  for  the  Lear  flight  of 
Figure  1  since  the  2-D  instrument  was  inoperative  during  that  time  period. 

The  good  agreement  between  the  exponential  number -density  fit  and  the  Dj^ 
assumption  with  the  experimental  results  give  an  added  degree  of  confidence  for 
the  determination  of  the  parameters  from  a  predicted,  full  distribution.  Since  Nq 
and  A  are  known  and  a  dm  has  been  established,  the  estimated  values  of  M  and  Z 
can  be  determined  by  finding  the  M  and  Z  contribution  from  the  extrapolated  part 
of  the  distribution  (dotted  line  in  Figures  1,  2  and  3)  and  adding  them  to  the  meas¬ 
ured  values  as 


ME  =  Mp  +  AM 


Zg  =  Zp  -f-  AZ  mm6  m  3  (5 

where  Mg  and  Zg  are  the  estimated  liquid  water  content  and  radar  reflectivity, 
Mp  and  Zp  are  the  values  derived  from  the  precipitation -probe  measurements 
and  AM  and  AZ  are  the  contributions  from  the  extrapolated  DT  to  DM  portion. 
The  equations  for  Mg  and  Zg  thus  become 


ME  * 


Md  dD  + 


DD  dD  g  m 


f  r 

Zg  =  J  Zg  dD  +  j  Zjj  dD  mm6  m’3 


Since  both  first  terms  are  known  quantities,  the  solving  of  these  equations  be¬ 
comes  a  two-part  process  where  Mp  and  Zp  are  used  in  Eq.  (22)  to  determine  the 
A  of  the  exponential  function  which,  in  turn,  is  used  in  either  Eq.  (12)  or  (17)  to 
find  Nq.  These  parameters  are  then  incorporated  into  the  second  term  of  Eqs. 

(26)  and  (27)  to  find  AM  and  AZ. 

Table  2  lists  the  Mp  and  Zp  that  were  derived  from  the  measurements  from 
the  precipitation  probe  for  the  three  sample  passes  compared  with  the  estimated 
Mg  and  Zg  from  the  extrapolated  spectra. 
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Table  2.  Comparison  of  Precipitation  Probe  Data  (Mp,  Zp)  with 
the  Estimated  Values  <ME,  Z£>  Determined  from  the  Extrapolated 
Distributions 


Mp 
g  m"3 

me 

g  m’3 

Factor 

Mg/Mp 

2P 

6  -3 

mm  m 

ZE 

mm6  m”3 

Factor 

ZE/Zp 

Lear 

0.3279 

0.3533 

1.077 

177.2 

268.3 

1.514 

MC-130E 
Pass  3 

0.3824 

0.4144 

1.084 

544.1 

845.  5 

1.554 

MC-130E 
Pass  5 

0.3062 

0.3137 

1.024 

336.4 

393.6 

_ i 

1.170 

In  these  three  samples  the  increase,  as  a  result  of  extrapolation,  was  only 
in  the  order  of  ^2  to  8  percent  in  M  but  —  17  to  55  percent  in  Z. 

Although  the  amount  of  M  and  Z  missed  because  of  instrument  truncation 
could  vary  considerably  in  different  situations,  it  is  clear  from  these  results  that 
it  is  the  Z-parameter  which  shows  the  largest  degradation.  In  the  case  of  a  slight  ly 
truncated  distribution,  the  measured  M  may  be  an  accepted  value,  but  if  the  sam¬ 
pled  hydrometeors  are  to  be  defined  in  a  relationship  of  aircraft -measured  M  and 
Z,  then  new  estimated  values  of  Z  have  to  be  incorporated. 

?.  REGRESSION  ANALYSIS  OF  MULTI  SAMPLE  TRUNCATED 
AIRCRAFT  DATA 

In  the  normal  analytical  procedure  used  at  AFGL,  the  multi-sample  spectral 
data,  acquired  by  the  Lear  and  MC-130E  along  particular  flight  paths  through 
cloud  and  storm  situations,  are  used  to  obtain  power -function  equations  that  math¬ 
ematically  define  hydrometeor  environments.  15,  16  ^one  by 

the  conventional  method  of  cross -plotting  the  logarithmic  values  of  M  and  Z  from 
the  multiple,  individual  samples.  The  least  square  analysis  of  this  field  of  data 
points  determines  the  line  of  best  fit  and  provides  the  coefficient  and  exponent  of 
the  power -function  regression  equation.  The  ultimate  goal  is  to  relate  M  and  Z  in 
a  regression  relationship  as 

M  =  a  Zb  (28) 


Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  41. 
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so  that  the  liquid  water  content  of  a  specific  hydrometeor  region  can  be  calculated 
from  a  given  radar  reflectivity. 


J  Two  methods  have  been  employed  in  past  work  to  acquire  M-Z  information; 

f  one  is  the  regression  of  the  dependent  variables  of  aircraft  derived  versus 

Z ^  and  the  other  is  the  regression  of  the  independent  variables  of  versus  the 
radar -measured  ZR.  Since  the  aircraft -spectral  M  and  Z  values  are  highly 
dependent  upon  the  a  and  0  in  the  i  to  D  conversion  (Eq.  (4))  and  are  also  suscep¬ 
tible  to  Instrument  truncation  effects,  the  first  method  has  an  inherent  uncertainty 
in  both  of  the  values  being  correlated.  Using  the  Z^  in  place  of  the  Z^  in  the 
correlation  removes  these  particular  uncertainties  from  one  parameter  of  the 
analysis  since  Z^  does  not  require  a  prior  knowledge  of  a  and  0  and,  of  course, 
is  not  limited  to  the  measuring  range  of  the  1-D  instruments.  Another  benefit  of 
using  ZR  is  the  ability  to  sample  the  defined  region  during  periods  when  the  air¬ 
craft  must  be  temporarily  absent. 

Studies  have  been  conducted  at  AFGL  on  the  variability  of  the  spectral  and 
Z^  caused  by  an  incomplete  knowledge  of  the  a  and  0  terms  of  the  I  to  D  conversion. 
For  example,  if  the  aircraft  were  to  sample  a  homogeneous  environment  of  pris¬ 
tine-crystal  type  and  the  a  and  0  of  Eq.  (4)  were  precisely  known,  then  the  meas¬ 
ured  radar  reflectivity  (Z^)  should  equal  or  very  closely  approximate  the  Z^ 

\  calculated  from  the  aircraft  data.  Ideal  conditions  such  as  this  are  seldom,  if 

E  ever,  found  in  storm  situations  and  usually  the  crystal  mixes  are  so  varied  that 

the  a  and  0  is  normally  determined  from  the  predominant  type.  The  result  is  that 
Z^  and  ZR  do  not  always  closely  agree;  thus  the  liquid  water  content  derived  from 
the  aircraft  measurements  (M^)  is  not  always  the  true  M  of  the  sample.  It  has 
been  demonstrated  that  by  using  a  spectral  parameter  k  in  the  regression  analysis, 
where  k  is  defined  as 

k  =  MA/SZA  g  mm*3  m'1-  5  (29) 


that  the  uncertainty  effect  of  the  I  to  D  relationship  of  Eq.  (4)  is  minimized  when 


deriving  the  M  of  a  sample. 


9,  16 


The  correlation  of  k  with  the  radar  values  Is  given  by 
bi 

k  '  al  ZR 

and  is  related  to  the  M  of  the  sample  by 
0.  5+b, 


M  ■  *1  ZR 


"1 


(30) 


(31) 
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The  aircraft  values  themselves  can  also  be  expressed  as 


k  = 


a2  ZA 


(32) 


and 


0.  5+b, 

^  s  a2 


(33) 


If  the  aircraft  data  used  in  Eqs.  (29)  through  (33)  are  taken  from  the  precipi¬ 
tation  probe  only,  as  in  the  following  discussion,  then  the  subscript  "A"  should  be 
replaced  with  a  "  F"  to  more  accurately  describe  the  process.  "A"  will  then  refer 
to  the  combined  values  from  both  the  cloud  and  precipitation  probes. 

When  the  aircraft  spectral  data  used  in  these  regression  analyses  stem  from 
inst rum entaily -truncated  samples,  then  it  is  apparent  that  the  data  must  be  cor¬ 
rected  or  compensated  to  take  into  account  the  effects  of  the  truncation.  The 
correction,  when  applied  to  the  two  regression  methods,  will  display  a  similar 
effect  on  both  of  the  k-Z  correlations  even  though  both  parameters  of  the  kp-Zp 
regression  will  be  altered  by  the  compensation  and  only  the  aircraft  k-values  in 
the  kp-Zp  regression  will  be  affected.  These  two  regression  cases  will  be  con¬ 
sidered  separately  in  the  following  discussion  and  the  methods  that  may  be  used 
to  correct  them  for  the  effects  of  truncation  will  be  described. 

An  example  of  the  kp-Zp  correlation  is  plotted  in  Figure  4  and  shows  the  53 
data  points  (4-second  averages)  of  the  precipitation-probe  measurements  from  the 
MC-130E,  pass  5  (Figure  3)  scattered  about  the  regression  line.  Any  kp-Zp 
relationship  of  this  nature  will  allow  a  corresponding  value  of  kp  to  be  calculated 
for  any  given  value  of  Zp. 

If  some  or  all  of  the  sampled  distributions  are  instrumentaliy  truncated,  as  in 
the  example,  the  analysis  (Eqs.  (32)  and  (33))  is  still  valid  for  that  specified 
limited  situation  but  does  not  reflect  the  condition  that  actually  exists.  One  method 
of  correcting  this  regression  is  to  use  the  Mp  and  Zp  of  each  individual  truncated 
distribution  in  Eq.  (22)  and  estimate  a  new  and  Z„  as  outlined  in  Section  2. 

New  regressions  can  then  be  calculated  based  on  the  estimated  parameters 


a3  ZE 


(34) 


and 


0.  5+b, 
ME  "  a3  ZE 

as  shown  in  Figure  5. 


(35) 


21 


Figure  4.  The  kp-Zp  Regression  of  the  53 
Truncated  Data  Points  (4 •Second  Averages) 
from  the  MC-130E.  Pass  5 


kE=.09lll  ZE-2944 


Figure  5.  The  kg-Zg  Regression  of  the  53 
Extrapolated  Data  Points  from  the  MC-130E, 
Pass  5 


1 


This  procedure,  of  course,  is  valid  only  when  the  multi-sample,  averaging 
time-interval  of  the  spectra  is  long  enough  and  the  data  are  of  sufficient  quality  to 
insure  that  the  individual  distributions  are  of  exponential  form.  As  a  check  of  this 
sample  meeting  these  criteria,  six  randomly -chosen  distributions  from  the  total 
of  53  are  plotted  in  Figure  6.  When  compared  with  the  pass  average  plot  (Figure  3), 
it  is  evident  that  the  much  shorter  averaging  time-interval  causes  an  increase  in 
the  variability  of  the  number  density  about  the  exponential  slope  although  each  of 
the  six  cases  can  still  be  reasonably  described  by  an  exponential  function. 

The  method  of  regressing  the  estimated  parameters  from  the  extrapolation  of 
each  individual  point  (Figure  5)  displays  a  more  negative  slope,  which  is  as 
expected,  when  one  reflects  upon  the  changes  that  result  from  extrapolation. 

It  may  be  helpful,  at  this  point,  to  review  the  changes  that  will  occur  when  the 
extrapolation  process  is  used  and  to  illustrate  in  a  simple  diagram  the  expected 
effects  of  these  corrections. 

Figure  7  describes  two  kp-Zp  regression  lines  (positive  and  negative)  which 
are  formed  from  five  truncated  distributions  (open  circles).  In  Section  2  it  was 
shown  that  the  number -density  plot  of  an  instrument  ally  truncated  distribution, 
when  extrapolated  to  a  new  D^A  limit,  will  have  a  much  larger  increase  in  Z  then 
in  M  since  M  varies  with  diameter  to  the  third  power  (Eq.  (1»  and  Z  varies  with 
diameter  to  the  six  power  (Eq.  (3)).  Also,  when  considering  a  distribution  with  a 
particular  exponential  slope,  the  larger  the  M  and  Z  in  the  truncated  distribution 
the  larger  the  change  will  be  upon  extrapolation.  Since  the  k-parapneter  is  defined 
as  M  divided  by  the  square  root  of  Z  (Eq.  (29»  then  it  follows  that  k  will  decrease 
as  a  result  of  extrapolation  and  will  also  show  larger  changes  in  the  more  densely 
populated  cases. 

The  X's  in  Figure  7  represent  the  estimated  k£  and 
from  correcting  each  truncated  distribution  and  show  the  decrease  in  k  along  with 
the  increase  in  Z  for  each  compensated  point.  The  resulting  changes  in  both  the 
positive  and  negative  regression  lines  are  similar,  in  that  they  both  exhibit  a 
more  negative  slope  because  of  the  corrections. 

The  method  of  extrapolating  each  individual  distribution  (Figure  5)  gives  satis¬ 
factory  results  but,  since  it  involves  so  many  extrapolations  (one  for  each  of  the 
multiple  samples),  it  becomes  a  long  complicated  process  and  it  is  obvious  that  a 
simpler  and  quicker  procedure  is  needed. 

A  regression  line  calculated  from  truncated  data,  such  as  in  Figure  4,  repre¬ 
sents  the  best  straight-line  relationship  of  kp  and  Zp  as  dictated  by  the  hydrome¬ 
teors  sampled  during  the  pass.  The  regression  line  of  Figure  5  is,  correspondingly, 
the  best  relationship  obtained  by  the  regression  of  the  extrapolated  data  points. 
Transformation  from  the  truncated  to  the  non-truncated  case  may,  therefore,  be 
considered  in  the  context  of  change  in  the  lines  of  best  fit  rather  than  changes  In 


Zg  values  that  will  result 
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Figure  7*  Illustration  of  the  Expected  Change  in  the  k„ Regression  as  a 
Result  of  Extrapolation 


individua  distributions.  This  results  In  a  simpler  procedure  since  a  line  can  be 
defined  by  just  two  points  and  allows  the  "best-fit"  values  of  kp  and  Zp  from  two 
points  on  the  truncated  regression  line  to  be  extrapolated  to  give  values  of  k^,  and 
Zg.  These  estimated  values,  when  connected  by  a  straight  line  on  a  log-log  plot, 
will  give  the  predicted  regression  line  of  the  extrapolated  data*  Thus,  for  any 
Zp  there  is  a  corresponding  average  Mp  which  can  be  found  using  Eq.  (33)  and 
then  applied  to  the  extrapolation  process  through  Eq.  (22)  as 


A  s  61.2 


mm 


-1 


(36) 


Mp  can  also  be  found  through  Eq.  (29)  as 


Mp  =  k  >Tzp  g  m“3 


(37) 


and  can  be  Incorporated  into  Eq.  (22)  as 


The  selection  of  two  Zp's  and  computation  of  the  corresponding  Mp's  defines 
two  such  conditions  (Zpl,  Mpl  and  Zp2,  Mp2)  and  allows  Z£1,  M£1  and  Z£2, 

ME2  to  k®  determined  by  the  extrapolation  procedure  of  Section  2.  With  these  two 
estimated  points,  the  predicted  regression  line  and  the  coefficients  and  exponents 
of  Eqs.  (34)  and  (35)  can  be  determined. 

This  simplified  method  was  used  with  the  regression  values  from  the  original 
truncated  data  from  the  MC-130E,  pass  5  where  Zpl  and  Zp2  were  given  the  values 
of  plus  and  minus  one  standard  deviation  from  the  log-average  Z  of  the  pass. 

Table  3  lists  the  result  of  this  predicated  regression  in  the  a  and  b  terms  of  the 
power- function  equation  along  with  the  results  from  the  analysis  of  the  original 
data  and  from  the  extrapolation  of  each  of  the  53  individual  distributions. 

A  comparison  of  the  results  show  that  the  two  methods  of  estimating  the  re¬ 
gression  of  a  full  distribution  agree  quite  well  and  indicate  that  the  simplified 
procedure  may  be  used  in  place  of  the  longer  process  with  a  very  slight,  acceptable 
deviation  in  the  a  and  b  terms. 

In  both  of  the  methods  just  described,  the  estimated  M£  may  be  determined  for 
any  given  Z£  by  using  Eq.  (35)  but  it  may  be  remembered  that  if  the  total  is 
required,  the  amount  measured  by  the  cloud  probe  (M^)  should  be  taken  into 
account. 
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Table  3,  Comparison  of  the  a  and  b  Terms  of  the 
Power-Function  Equation  from  the  Original  Data, 
Extrapolation  of  Each  Separate  Distribution  and 
Predicted  Regression 


a 

b 

MC-130E,  pass  5  original  data 

0.  08203 

-0.2722 

Each  of  the  53  separate  dis¬ 
tributions  extrapolated 

0.09111 

-0.  2944 

Predicted  regression 
(Simplified  Method) 

0.09206 

-0.  2931 

Discussion  of  the  cloud  probe  data  and  its  combination  with  that  from  the 
precipitation  probe  is  contained  in  Appendix  A, 

When  a  particular  hydrometeor  environment  is  measured  with  both  radar  and 
aircraft  and  the  data  from  the  two  sources  are  correlated,  then  Eq.  (30)  is  used 
to  define  the  k  that  is  associated  with  ZR.  In  situations  where  the  aircraft  meas¬ 
urements  are  instrument  ally  truncated,  it  is  the  k-parameter  only  that  will  be 
altered  by  the  application  of  the  extrapolation  procedure  since  the  ZD  measure - 

rt 

ments  are  not  subject  to  the  truncation  limits  of  the  1-D  instrument.  Therefore, 
the  values  of  ZR  are  the  actual  reflectivity  factors  from  the  full  distributions.  A 
kp-ZR  regression  line  produced  through  the  correlation  of  truncated  aircraft  data 
with  radar  measured  ZR's  will  exhibit  a  different  slope  than  the  kp-Zp  regression 
line  but  should  change  (assuming  that  ZE  is  equal  to  ZR),  upon  extrapolation,  to 
match  the  slope  of  the  k£-ZE  regression.  (In  actual  practice,  Z£  should  not  be 
expected  to  exactly  equal  ZR  since  both* the  aircraft  and  radar  measurements  con¬ 
tain  inherent  instrumentation  and  data-acquisition  errors. ) 

Once  again,  an  illustration  may  be  helpful  in  visualizing  this  predicted  change. 
In  Figure  8  the  diagrams  of  Figure  7  have  been  altered  to  include  the  kp-ZR  trun¬ 
cated  data  points  (closed  circles)  and  their  direction  of  change  upon  extrapolation 
for  regression  lines  of  positive  and  negative  slopes.  In  both  cases,  the  kp-ZR 
regression  lines  are  closer  to  a  zero  slope  than  their  kp-Zp  counterparts  but  the 
most  noticeable  difference  occurs  in  the  amount  of  change,  with  the  k„-Zn  line 

r  x\ 

varying  more  in  the  negative  case  and  the  kp-Zp  line  showing  the  greatest  change 
in  the  positive  case. 

In  experimental  situations  where  aircraft  and  radar  data  are  available  and  the 

aircraft  values  are  from  truncated  measurements,  it  is  required  to  compensate 

kp.  In  order  to  Incorporate  the  extrapolation  process  for  the  correction  of  kp,  it 

is  necessary  to  determine  the  Mp  and  Zp  for  any  given  ZR.  This  is  easily 

accomplished  by  forming  a  new  regression  of  Z_,  and  Z0  as 

P  n 
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(39) 


r 


z 


p 


which,  when  substituted  in  Eq.  (29),  gives 


This  allows  Eq.  (22)  to  be  expressed  in  terms  of  kp  and  ZR  as 


(40) 


(41) 


Using  the  same  technique  as  in  the  simplified  method  of  the  kp-Zp  case,  two 
selections  of  Z  may  be  chosen,  this  time  Z^^  and  Zpg»  alM*  through  Eq.  (41)  and 
the  extrapolation  process  described  in  Section  2,  two  new  k^s  can  be  evaluated. 
These  k^'s,  with  the  selected  ZR‘s,  will  give  the  estimated  relationships  of  the 
non -truncated  spectra  with  the  radar  data  as 


k 


E 


(42) 


and 


0.  5+bj. 

ME  =  a5  ZR  5  (43) 

The  reader  is  again  referred  to  Appendix  A  for  a  discussion  on  the  combina¬ 
tion  of  the  cloud  and  precipitation  data. 


4.  TEST  OF  THE  REGRESSION-PREDICTION  METHODS 

In  Section  2,  it  was  demonstrated  than  an  instrument  ally -truncated  distribution 
can  be  defined  by  an  exponential  function  when  converted  into  terms  of  equivalent 
melted  diameters.  Therefore,  to  test  the  regress  ion -predict  ion  methods  outlined 
in  Section  3,  exponential  distributions  (non -instrument ally  truncated)  may  be  de¬ 
fined  as  standards  and  then  truncated,  by  reduction  of  the  limits,  to  be  made  to 
represent  instrument  ally  truncated  aircraft  measurements.  The  values  calculated 
from  these  truncated  situations  can  then  be  processed  through  the  previously 
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described  regress  ion -predict  ion  procedures  and  the  results  compared  with  those 
of  the  full  distributions. 

Four  exponential  distributions  were  chosen  so  that  they  would  have  differing 
slopes  and  that  their  k  and  Z  points  would  scatter  about,  and  not  fall  directly  on 
the  regression  line.  In  all  of  the  four  cases,  the  minimum  diameter  was  0. 1  mm 
and  the  maximum  was  2.  6  mm.  The  parameters  of  these  spectra  are  listed  in 
Table  4  and  their  number -density  plots  are  shown  in  Figure  9. 


Table  4.  Parameters  of  the  Four  Distributions  Used  to  Test 
the  Regression-Prediction  Methods 


When  the  k's  and  Z's  from  Table  4  are  subjected  to  regression  analysis  as  in 

-0  00379 

Figure  10,  the  power -function  equation  is  k  =  0.010085  Z  *  and  can  be 

designated  as  the  standard  equation  to  which  any  predicted  regression  can  be  com¬ 
pared.  Then  each  of  the  test  spectrums  can  be  truncated  at  any  diameter  (DT>, 
where  DT  is  always  less  than  DM,  and  the  k  and  Z  calculated  under  this  condition 
will  be  synonymous  with  that  derived  from  truncated  aircraft  data. 
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Figure  10.  The  k-ZR  Regression  from  the  Four  Test 
Distributions  ** 


This  was  done  for  each  of  the  samples  for  the  DT  limits  of  1.  0,  1.2,  1.4,  1.6, 
1.8,  2.  0,  2.  2  and  2.  4  mm  and  the  resulting  k's  and  Z's  are  listed  in  Table  5. 

At  this  point,  the  simulated  aircraft  data  has  been  defined  and  consideration 
will  now  be  given  to  the  application  of  these  data  to  the  previously  outlined  methods 
of  predicting  the  regression  equation  from  the  truncated  distributions. 

The  first  method  discussed  in  Section  3  was  that  of  aircraft  values  alone.  In 
this  case,  the  k  and  Z^  of  each  distribution,  that  are  the  result  of  a  particular 
truncation  limit,  were  used  in  a  regression  analysis  to  give  the  relationships  of 
Eqs.  (32)  and  (33).  Two  Z^'s  were  then  designated  as  previously  described  (plus 
and  minus  one  standard  deviation  from  the  log  average)  and  these  Z^’s  were  then 
applied  to  Eq.  (36)  and  the  extrapolation  process  of  Section  2.  The  resulting  Mg's 
and  Zg's ,  one  pair  for  each  of  the  two  defined  Z^'s,  were  then  used  to  give  the 
predicted  regression  lines  of  the  extrapolated  data  and  the  coefficients  and  expo¬ 
nents  of  the  power -function  relationship  of  Eqs.  (34)  and  (35). 
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Table  5.  The  Simulated  k.  fs  and  Z.'s  of  the  Test  Distributions  Under 
Truncated  Conditions 


Dt  mm 

Dist.  1 

Dist.  2 

Diet.  3 

Dist.  4 

1  A 

kg  mm"3  m“l*  5 

Z  mmD 

0.014979 

0.016258 

0.010745 

0.011830 

1.  u 

103.481 

71.058 

25.  032 

10.232 

1.2 

k 

0.013534 

0.014625 

0.009718 

0.011165 

Z 

190.266 

107.  564 

35.  698 

11.910 

1.4 

k 

0.012353 

0.013509 

0. 009060 

0.010865 

Z 

292.347 

139.  169 

43.990 

12.  707 

1.6 

k 

0.011427 

0.012781 

0. 008660 

0.010742 

Z 

396.473 

162.  936 

49.571 

13.034 

1.8 

k 

0.010723 

0.012326 

0. 008430 

0. 010698 

Z 

491.889 

178.930 

52. 942 

13. 154 

El 

k 

0.  010200 

0.012054 

0. 008305 

0.010681 

Z 

572.316 

188. 837 

54.814 

13. 195 

9  9 

k 

0. 009822 

0.011899 

0. 008239 

0.010676 

6.  C 

Z 

635.  721 

194.575 

55.  786 

13.208 

wgm 

k 

0. 009555 

0.011813 

0. 008208 

0. 010675 

mam 

Z 

683.042 

197.719 

56.263 

13.211 

f 


f 

i 

f 

j- 

F 


This  procedure  was  followed  in  each  of  the  truncated  cases  and  Table  6  lists 
the  coefficients  and  exponents  of  the  power -function  equations  for  each  DT  limit 
along  with  those  for  the  standard  equation  from  the  full  distributions. 

The  regression  lines  corresponding  to  each  Dj,  are  plotted  in  Figure  11  and 
show  the  negative  change  in  slope  as  approaches  the  Dj^  limit,  (dotted  line) 
Table  7  gives  the  coefficients  and  exponents  of  the  predicted  regression  equations 


and  compares  them  with  those  from  the  standard  equation. 

It  is  apparent  that  there  are  only  very  slight  differences  between  the  predicted 
equations  and  the  standard,  but  since  M  is  the  final  product,  the  evaluation  can  not 
be  considered  complete  until  the  M  values  from  the  predicted  and  standard  equa¬ 


tions  are  compared.  This  was  done  for  designated  Z!s  of  10,  100  and  1000  and  the 


results  are  listed  in  Table  8. 


The  second  prediction  method  described  in  Section  3  was  the  combination  of 
aircraft  and  radar  data.  To  test  this  method,  the  radar  is  given  the  Z  values  of  the 
non -truncated  distributions  from  Table  4.  The  truncated  aircraft  values  are,  once 
again,  those  listed  in  Table  5. 
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Table  6.  The  Coefficients  and  Exponents  from  the  Truncated 
and  Non -Truncated  Cases 


DT  mm 

k  =  a2Z 

^  —  M  =  a2 

0.  5+b0 

ZA 

a2 

b2 

0.5 +  b2 

1.0 

0.007767 

0.14816 

0.64816 

1.2 

0.007891 

0. 10714 

0.  60714 

1.4 

0. 008261 

0. 07440 

0. 57440 

1.6 

0.008692 

0.04941 

0.54941 

1.8 

0.009100 

0. 03083 

0.  53083 

2.0 

0.  009445 

0.01734 

0.51734 

2.2 

0.009721 

0.00769 

0. 50769 

2.4 

0. 009932 

0.  00087 

0. 50087 

%  2.6 

0.010085 

-0.  00379 

0.49621 

Table  7.  Predicted  Regressions  and  Their  Comparison  with  the  Standard 
(Simulated  Aircraft) 


Dt  mm 

t*T 

ii 

p» 

09 

N 

b3 

E  ME  =  a2 

0.  5+b, 
lZE 

Factor  =  Predlcted 
*  actor  Standard 

a3 

b3 

0.  5  +  b3 

a 

0.  5  +  b 

1.0 

0.010067 

•0. 00334 

0. 49666 

0. 9982 

1.0009 

1.2 

0, 010008 

-0.  00344a 

0.49656 

0.9924 

1.0007 

1.4 

0.010211 

-0. 00680 

0.49320 

1.0125 

0.  9939 

1.6 

0.010150 

-0. 00632 

0. 49377 

1.0064 

0.  9951 

1,8 

0. 010086 

-0. 00405 

0.49595 

1.0000 

0.9995 

2.0 

0.010051 

-0. 00320 

0.49680 

0.9966 

1.0011 

2.2 

0.010046 

-0. 00302 

0. 49698 

0.9961 

1.0015 

2.4 

0.010061 

-0. 00331 

0.  49669 

0.9976 

1.0010 

DM2.6 

0. 010085 

-0.00379 

0. 49621 

i 

1 _ 

33 


Z  mm6 


Figure  11.  The  k-Z^  Regressions  from  the  Four 
Test  Distributions  at  Various  Truncation  Limits  (D^) 


Table  8.  Comparison  of  M  from  the  Predicted  and  Standard  Equations 
(Simulated  Aircraft).  Factor  (F)  -  predicted /standard 


Dt  mm 

Z  =  10  mm  3  m  3 

Z  =  100  mm6  m  3 

Z  =  1000  mm6  m  3 

_3 

M  g  m 

F 

..  "3 

M  g  m 

F 

M  g  m"3 

F 

1.0 

0.031591 

0. 9993 

0.099133 

1.0003 

0.311086 

1.0013 

1.2 

0.031398 

0. 9932 

0.098507 

0.9940 

0.309049 

0. 9948 

1.4 

0.031788 

1.0055 

0.098962 

0.9986 

0.308083 

0.9917 

1.  6 

0.031640 

1.0008 

0. 098629 

0.9952 

0.307451 

0. 9896 

1.  8 

0.031599 

0. 9995 

0.098996 

0.9989 

0.310148 

0. 9983 

2.0 

0,  031550 

0.  9980 

0.  099040 

0. 9993 

0.310892 

1.0007 

2.2 

0.031548 

0.9979 

0.099072 

0.9997 

0.311124 

1.0014 

2.4 

0.031574 

0. 9987 

0.099088 

0. 9998 

0.310965 

1.0009 

dm2*6 

0.031614 

— 

0.099105 

— 

0.310675 

-- 

In  this  method,  both  the  k  and  Z^  of  each  distribution,  that  are  the  result  of 
a  specific  truncation  limit,  are  related  to  the  assigned  ZR  to  give  the  k-ZR 
(Eq.  (30))  and  Z^-ZR  (Eq.  (39))  regressions.  Two  ZR's  are  then  designated 
(once  again,  plus  and  minus  one  standard  deviation  from  the  log  average)  and  the 
k’s  associated  with  the  two  ZR  values  are  determined  for  the  specified  DT  limit. 
Substituting  the  two  ZR's  and  their  associated  k's  in  Eq.  (41)  then  allows  two 
estimated  (or  k^)  and  Z£  pairs  to  be  calculated  through  the  extrapolation 
process  of  Section  2.  The  plotting  of  the  logarithmic  values  of  these  estimated 
points  will  give  the  predicted  regression  line. 

This  procedure  was  followed  for  each  of  the  truncated  situations  and  the 
coefficient  and  exponents  of  the  k-ZR  and  Z^-ZR  relationships  are  listed  in 
Table  9. 


Table  9.  The  Coefficients  and  Exponents  of  the  k-ZR  and 
Z^-ZR  Regressions  K 


Dt  mm 


7? 

ii 

zbi 

ZR 

b4 

ZA  =  a4  ZR 

al 

bl 

a4 

b4 

0. 009023 

0.08338 

2.27332 

0.  60392 

0. 008646 

0. 07284 

2.02020 

0.71226 

0.008682 

0. 05745 

1.  72348 

0.  79862 

0.008911 

0.04152 

1.47858 

0.  86462 

0.009209 

0.02734 

1.29992 

0.91327 

0. 009497 

0.01586 

1.  17622 

0.94807 

0.009744 

0. 00711 

1.09254 

0.97235 

0.009941 

0. 00069 

1.03664 

0. 98893 

The  predicted  regressions  are  given  in  Table  10  where  they  are  compared  with 
the  standard  equation. 

Once  again,  the  more  significant  comparison  of  M's  from  the  predicted  and 
standard  equations  are  given  in  Table  11. 

It  is  evident,  that  the  factors  in  the  radar-aircraft  case  have  a  slightly  larger 
variance  than  that  of  the  aircraft -only  method,  although  both  are  able  to  estimate 
the  standard  regression  equation  to  within  a  very  close  tolerance. 


Table  10.  Predicted  Regressions  and  Their  Comparison  with  the 
Standard  (Simulated  Aircraft  and  Radar) 


D,p  mm 

mm 

HI 

a5 

b5 

HHB9 

a 

| 

1.0 

0. 010462 

-0. 01184 

0.48816 

1.0374 

0.  9838 

1.2 

0. 010420 

-0.01112 

0. 48888 

1.0332 

0.  9852 

1.4 

0.010321 

-0.00913 

0. 49087 

1.0234 

0.  9892 

1.6 

0.  010208 

-0. 00676 

0.49324 

1.0122 

0. 9940 

1.8 

0.010127 

-0. 00494 

0.49506 

1.0042 

0.  9977 

2.0 

0. 010079 

-0. 00380 

0. 49620 

1.0000 

2.2 

0. 010063 

-0.  00339 

0.49661 

■ 

1.0008 

2.4 

0. 010068 

-0. 00344 

0.49656 

■ 

1.0007 

°M  2#6 

0.  010085 

-0.00379 

0.49621 

| 

-- 

Table  11.  Comparison  of  M  from  the  Predicted  and  Standard  Equations 
(Simulated  Aircraft  and  Radar).  Factory  (F)  =  predicted /standard 


Z  -  10  mm**  m”3 

Z  -  100  mr 

6  -3 

n  m 

mol 

Dt  mm 

M  g  m~3 

F 

M  g  m“3 

F 

M  g  m  3 

F 

1.0 

0.032194 

1.0183 

0. 099068 

0. 9996 

0.304856 

0.  9813 

1.2 

0.032118 

1.0159 

0. 098998 

0.9989 

0.305146 

0.  9822 

1.4 

0.031959 

1.0109 

0. 098960 

0.  9985 

0.306430 

0.  9863 

1.6 

0. 031782 

1.0053 

0. 098951 

0. 9984 

0.308078 

0.  9916 

1.8 

0. 031662 

1.0015 

0.  098992 

0. 9989 

0.309500 

0. 9962 

2.0 

0.031595 

0.9994 

0.  099042 

0.9994 

0.310468 

0. 9993 

2.2 

0.031575 

0. 9988 

0.099071 

0.9997 

0.310855 

1.0006 

2.4 

0. 031866 

0.9991 

0. 099098 

0. 9999 

0.310902 

1.0007 

dm2-6 

0.031614 

-- 

0.099105 

— 

0.310675 

-- 
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5.  REAL-TIME  ESTIMATIONS  FROM  TRUNCATED  DISTRIBUTION 
DATA 


Real-time  analysis  of  the  1-D  data  is  very  important  to  the  meteorologist  on 
board  the  sampling  aircraft  and  also,  through  a  radio  link,  to  the  ground-based 
mission  control.  Knowledge  of  liquid  water  content  can  be  particularly  critical  in 
some  situations  especially  in  the  lower  altitudes  of  a  storm  where  snow  crystals 
are  large  and  truncation  effects  are  most  likely  to  occur. 

Both  the  Lear  and  the  MC-130E  have  on-board  computers  which  give  the  field 
personnel  real-time  parameters  of  the  hydrometeor  environment  in  which  they  are 
sampling.  However,  when  the  aircraft  passes  through  a  region  containing  large- 
size  particles  the  calculated  results  reflect  the  conditions  of  truncated  spectra 
and  not  the  true  situation  of  the  full  actual  distributions.  It  is  at  this  point  that 
quick,  yet  reasonably  accurate  estimates  of  the  true  M  and  Z  are  needed. 

It  is  a  relatively  simple  matter  to  instruct  the  computers  to  recognize  a  prob¬ 
able  truncated  situation  since  this  possibility  could  exist  whenever  the  largest 
diameter-size  channel  of  the  precipitation  probe  has  recorded  counts.  But,  it  is 
not  feasible  to  have  these  computers  extrapolate  distributions  because  of  their 
limited  capacity. 

Therefore,  equations  have  been  derived  to  give  an  approximation  of  the 
and  Zg  that  would  be  obtained  through  the  extrapolation  of  a  truncated  distribution 
using  as  input  the  real-time  calculations  of  Mp  and  Zp.  These  equations  are  based 
on  exponential  distributions  that  have  upper  and  lower  diameter  limits  correspond¬ 
ing  to  those  of  the  precipitation  probes  on  the  aircraft.  They  were  developed  using 
a  process  of  plotting  variations  of  Mp/Zp  vs  M^/Z^  until  a  very  near  straight- 
line  relationship  existed  and  then  were  converted  into  equation  form. 

A  large-snow  crystal  type  was  assumed  where  a  =  0.  4  and  0  =  0.  782  in 
Eq.  (4)  and  the  extrapolation  limit  was  set  at  E>MA  =  9. 

The  equations  for  large  snow  conditions  using. the  Lear  probe  are 


Me 


M, 


2303 

e 


Mp  X  10 


•9.  5155 


0.  02 


and 


ZE  '  ZP 


394.  3 
e 


Mp  X  10 


a  v-6.  1093 

-)  ...» 


(44) 


(45) 


I 


i 
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The  cor  responding  equations  to  be  used  with  the  probe  on  the  MC-130E  are 


(46) 


(47) 


Both  sets  of  equations  are  of  the  same  form  with  the  exponents  reflecting  the 
differences  in  the  upper  truncation  limits  of  the  two  different  1-D  precipitation 
probes. 

It  should  be  emphasized,  that  these  equations  have  been  developed  to  give  a 
quick  method  of  estimating  the  M  and  Z  irom  a  single  distribution  (that  is,  aircraft 
data  averaged  over  any  length  of  time  into  a  single  distribution)  and  should  not  be 
confused  with  the  more  precise  post-mission,  multi-sample  regression  analysis. 

The  Mp's  and  Zp's  from  the  sample  distributions  presented  in  Section  2  were 
substituted  in  these  equations  and  the  results  are  listed  in  Table  12. 


Table  12.  Comparison  of  and  Z^  from  Extrapolation  and  Real-Time 
Equations 


Real 

-Time 

Ext  rapolated 

Factor 

/  Real-Time  \ 
^Extrapolated  J 

Me 

ZE 

me 

ZE 

g  m*3 

mm6  m3 

g  rn"3 

mm6  m~3 

for  Me 

for  ZE 

Lear 

0.3  509 

268.  5 

0.3533 

268.3 

0.993 

1.001 

MC-130E 
Pass  3 

0.4162 

844.2 

0.4144 

845.  5 

1.004 

0.998 

MC-130E 
Pass  5 

0.3159 

393.4 

0.3137 

393.6 

1.007 

1.000 

Further  comparisons  of  and  Z£  wer^  made  between  the  results  of  the 
real-time  calculations  and  the  extrapolation  process  using  a  variety  of  Mp's  and 
Zp*s.  It  was  found,  through  this  testing,  that  if  the  or  ZE  were  less  than 
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twice  the  value  of  Mp  or  Zp,  then  the  estimated  results  from  the  two  methods 
differed  less  than  1  percent.  In  other  words,  if  the  estimated  M  were  within  a 
factor  of  2  from  Mp  then  the  results  calculated  by  either  method  fall  within 
1  percent  of  each  other.  The  same  rules  apply  for  Zp  and  Z£. 

When  the  tests  were  compared  at  a  factor  of  3  increase  in  the  estimated 
values,  the  differences  were  within  5  percent . 

As  a  final  word,  it  must  be  remembered  that  the  subscripted  "P"  and  "E" 
values  refer  to  the  measured  and  estimated  data  of  the  precipitation  probe  and 
that  the  and  Z^  from  the  cloud  probe  must  be  added  to  the  estimated  values 
to  obtain  a  total  and  Z^. 

6.  CONCLUSIONS 

The  analytical  processes  presented  in  this  report  are  the  results  of  a  study 
conducted  on  truncated  measurements  from  aircraft  (PMS  1-D)  instruments.  A 
method  of  extrapolating  these  truncated  data  is  demonstrated  and  the  procedure 
whereby  the  extrapolated  data  are  conveniently  used  in  a  regression  analysis  is 
described. 

Testing  of  these  processes  have  indicated  that  this  extrapolation  method  can 
be  considered  valid  if  the  assumptions  of  an  exponential  number -density  distribution 
and  a  DM  upper  limit,  as  determined  by  the  equation  D^A  =  C,  are  accepted  as  a 
reasonably  accurate  description  of  reality  for  the  majority  of  cases.  However,  it 
must  be  emphasized  that  not  all  samples  will  have  distributions  that  are  able  to  be 
described  by  an  exponential  function.  An  occasional  pass  will  have  an  excess  or 
a  deficiency  of  specific  size  particles  which  result  in  number-density  distributions 
that  depart  from  exponential.  Also  the  D^A  constant  of  9  used  in  this  report  is 
specific  for  the  ice  hydrometeors  that  were  measured  in  a  particular  storm  and 
can  vary  with  crystal  shape  and  density.  Past  AFGL  measurements  have  indicated 
aDMA  constant  within  the  range  of  9  to  12. 

Instruments  of  the  future  may  prove  or  disprove  these  assumptions  and  may, 
ideally,  negate  the  need  of  extrapolating,  but  further  study  is  needed  to  related 
maximum  particle  crystal  size,  through  i  to  D  conversion,  to  a  DM  limit  until 
such  instruments  are  developed  and  put  into  operation,  A  better  knowledge  of  the 
dma  constant  and  the  or  and  0  of  the  I  to  D  equation  should  enable  the  exponent ial - 
distribution  equations  to  be  fine-tuned  to  more  accurately  describe  the  various 
snow  types  as  to  their  density  and  size. 

Although  the  described  extrapolation  method  deals  exclusively  with  ice  hydro - 
meteors,  it  can  also  be  applied  in  cases  of  truncated  distributions  in  rain.  The 
Lear  1-D,  for  example,  has  a  maximum  measuring  size  of  3.  105  mm  for  water 
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drops  and  can  be  expected  to  encounter  situations  having  hydrometeors  of  some 
larger  size.  In  cases  involving  liquid  water  there  is  no  i  to  D  considerations  and 
the  only  change  in  the  analytical  procedure  would  be  in  the  D^A  constant.  For 
rain,  this  constant  is  given  the  value  of  7.  5  up  to  an  absolute  maximum  of  5  mm 
which  is  considered  to  be  the  breakup  size  of  a  water  drop. 

It  is  recommended  that  the  methods  outlined  in  this  report  be  incorporated  in 
the  analyses  of  PMS  aircraft  data  and  that  the  real-time  equations  be  programmed 
into  aircraft  computers  so  that  full  utilization  may  be  made  of  the  instrument  ally 
truncated  data. 
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Appendix  A 

Th»  Combining  of  the  Cloud  and  Precipitation  Probe  Data 


The  extrapolation  process  described  in  this  report  is  concerned  only  with  the 
number -density  spectra  from  truncated  distributions  as  measured  with  the  1-D 
precipitation  probe.  Any  estimated  values  determined  through  the  use  of  this 
method  may  be  supplemented  with  the  data  from  the  cloud  probe  in  order  to  define 
the  parameters  of  a  total  spectrum  measured  over  the  range  of  both  instruments. 

Unlike  the  precipitation  probes,  the  cloud  1-Dfs  are  identical  on  both  aircraft 
with  each  being  capable  of  measuring  snow  crystals  from  0.015  to  0.360  mm. 
Since  the  smallest  measurable  sizes  of  the  precipitation  probes  are  0. 151  mm  for 
the  Lear  and  0.227  mm  for  the  MC-130E,  it  is  evident  that  an  overlap  occurs  in 
the  measurements  when  using  data  from  both  the  cloud  and  precipitation  probes. 
Fortunately,  the  redundant  spectra  can  be  conveniently  eliminated  by  considering 
just  the  cloud  instrument  data  for  all  sizes  up  to  that  which  corresponds  to  the 
lower  limit  of  the  precipitation  probe.  For  example,  on  the  Lear,  the  upper  limit 
of  channel  6  of  the  cloud  probe  is  0.  153  mm  and  when  compared  to  the  lower  limit 
of  0. 151  mm  of  the  precipitation  probe,  shows  a  slight  overlap  of  0.  002  mm. 
Channel  10  of  the  cloud  probe  has  an  upper  limit  of  0.  222  mm  and  when  compared 
with  the  lower  limit  of  0.227  mm  of  the  MC-130E  precipitation  probe  shows  that  a 
0.  005  mm  segment  is  missed.  These  errors  are  minor  and  can  be  safely  ignored. 

In  the  following  discussion,  all  references  to  cloud  probe  data  assume  that 
these  data  have  been  limited  in  accordance  with  the  above  rules. 
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Table  1  lists  the  M  and  Z  contributions  from  the  cloud  and  precipitation 
probes  for  the  three  sample  truncated  distributions.  It  is  evident  from  these  data 
that  when  extrapolation  of  the  data  is  required,  the  combined  cloud  and  precipi¬ 
tation  measurements  need  only  be  considered  for  the  M  parameter  since  the  cloud 
Z-contributions  are  0.01  percent  or  less  of  the  total.  In  some  cases,  which  de¬ 
pend  entirely  upon  the  proposed  use  of  the  data,  even  the  —  6  to  7  percent  M  con¬ 
tribution  may  be  of  such  little  importance  that  it  too,  may  be  ignored. 

The  combining  of  the  data  from  the  two  size  ranges  is  simple  and  straight¬ 
forward  when  considering  a  single  distribution.  In  this  instance,  the  is  added 
directly  to  th£  Mp  or  ME  to  give  the  total  as 

MA  =  MC  +  M(PorE)  <Al> 


When  two  or  more  distributions  are  used  to  form  a  kp-Zp  regression  any 
averaged  value  of  Zp  will  not  only  have  an  average  Mp  but  also  an  average 
associated  with  it.  This  can  be  determined  by  forming  a  relationship  between 
the  Mc  and  Zp  from  the  experimental  data  as 

Mc=xZyp  (A2) 

When  the  Zp  is  replaced  by  Z£  through  the  extrapolation  process  then  a  new 
relationship  can  be  formed. 


*1 

MC  =  X1  ZE 


(A3) 


The  determined  in  Eq.  (A3)  can  then  be  substituted  in  Eq.  (Al)  to  find  M^. 
Radar  Zrs  can  be  directly  associated  with  as 

Mc  =  x2  Zp2  <A4i 

where  the  is  then  used  in  Eq.  (Al). 

An  alternative,  that  may  be  used  to  avoid  these  regressions,  is  to  use  the 
pass-averaged  in  Eq.  (Al).  Since  the  percentage  of  the  averaged  is  small, 
when  compared  with  the  pass -averaged  total  the  discrepancies  caused  by  the 
direct  addition  of  the  mean  with  Mp  are  insignificant. 

These  two  methods  of  combining  the  cloud  and  precipitation  data  were  explored 
using  the  experimental  results  of  the  MC-130E,  pass  5  (Table  1,  Figures  3  and  4). 
This  pass  has  the  largest  contribution  of  (7.4  percent)  of  ail  three  samples 
listed  in  Table  1. 
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The  and  Z_  values  for  the  53  distributions  were  analyzed. and  gave  the 
i-  r  -0  1416 

relationship  of  Mc  =  0.  05569  Zp  *  .  The  maximum  Zp  of  the  pass  was  638.  6 

and  the  minimum  was  174.  2  and  when  applied  to  the  M^,-Zp  equation,  give  M^*s 
of  0.  0223  for  the  maximum.  Zp  and  0. 0268  for  the  minimum. 

The  equation  from  the  precipitation  probe  data  (Figure  4)  was  kp  =  0.  08203 
Z^0. 2722  when  converted  to  the  M-Z  form  of  Eq.  (33),  gives  an  Mp  of  0. 3573 
for  the  maximum  Zp  and  0.  2658  for  the  minimum. 

Combining  the  Mp's  with  the  Mc*s  from  the  regression  method  yields 
of  0. 3796  for  the  Zp  maximum  case  and  0.  2926  for  the  minimum.  The  pass- 
averaged  method  gives  0.3827  and  0.  2912  respectively  and  differs  from  the 
regression  values  by  less  than  1  percent. 

These  small  variations  justify  the  use  of  the  pass-averaged  values  of  Mc  in 
the  computation  of  total,  estimated  M  when  extrapolation  of  the  spectra  above  the 
upper  limit  of  the  1-D  precipitation  probe  is  required. 
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List  of  Symbols 


a  power -function  coefficient 

b  power-function  exponent 

C  constant 

d  minimum  drop  or  equivalent  melted  diameter 

drop  or  equivalent  melted  diameter  of  classified  data 
maximum  drop  or  equivalent  melted  diameter 
drop  or  equivalent  melted  diameter  at  truncation  limit 
k  aircraft  spectral  parameter 

kE  estimated  spectral  parameter 

I  measured  physical  size  of  ice  hydrometeors 

M  liquid  or  ice  water  content 

liquid  or  ice  water  content  from  aircraft  data 
liquid  or  ice  water  content  derived  from  cloud  probe 
Mp  liquid  or  ice  water  content  at  drop  or  equivalent  melted  diameter  D 

estimated  liquid  or  ice  water  content  from  an  extrapolated  distribution 
Mp  liquid  or  ice  water  content  derived  from  precipitation  probe 


AM  liquid  or  ice  water  content  from  the  extrapolated  part  of  the  distribution 

N  number  concentration  of  the  hydrometeor  drops  or  particles  per 
per  mm  bandwidth 

number  concentration  at  diameter  D 

number  concentration  for  classified  data 

Nq  the  "zero  intercept"  of  a  distribution  function  of  exponential  type 

Nt  total  number  of  the  hydrometeor  drops  of  particles  of  a  given  population, 
1  per  m^ 

r^  truncation  ratio  of  liquid  water  content 

r^  truncation  ratio  of  number  concentration 

r^  truncation  ratio  of  radar -reflectivity  factor 

Z  radar -reflectivity  factor 

radar -reflectivity  factor  computer  from  aircraft  data 
Z^,  radar-reflectivity  factor  from  cloud  probe 

ZD  radar -reflectivity  factor  at  diameter  D 

Zg  estimated  radar -reflectivity  factor  from  an  extrapolated  distribution 
Zp  radar -reflectivity  factor  from  precipitation  probe 
Zp  radar -reflectivity  factor  as  measured  by  radar 

AZ  radar -reflectivity  factor  from  the  extrapolated  part  of  the  distribution 
a  coefficient  of  the  i  to  D  equation 
0  exponent  of  the  I  to  D  equation 
r(n)  gamma  function  of  the  particular  number  "n" 

A  exponential  "slope  factor"  in  the  distribution  function  for  the  number 
concentration  of  the  drops  or  particles 

pw  density  of  liquid  water 
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